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(g) Discrete fourier transform circuit 

@ In a transmultiplexor, which is an interface between a 
number of PCM highways and an FDM group, the conver- 
sion between PCM and FDM involves first converting the 
PCM words Into their linear form and filtering them, apply- 
ing th m to a circuit which performs a discrete Fourier 
transf rm (DFT) on them followed by a second stage of 
filtering. An efficient implementation of a DFT is described 
which reduces the number of multiplications required by 
making use of certain symmetry properties in a DFT matrix. 

The whole can be used -backwards**, as it were to give 
an inverse Fourier transform as needed for the FDM-PCM 
conversion. 
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DISCRETE FOURIER TRANSFORM CIRCUIT 

This invention relates to an electrical circuit 
arrangement for implementing a discrete Fourier transform 
(DFT) , and also for implementing an inverse Fourier 
transform (IFT) • 

Such circuits are usable in digital signal 
processing applications and in one application a circuit 
comprising digital filters and a DFT implementation is 
used to provide the equivalent of -a set of filters used 
to assemble a frequency division multiplex (FDM) group 
from a number of PCM channels after the latter have been 
converted into a linear form. After the signals have 
been passed through the digital filtering arrangement 
thus formed they are passed to the FDM output via digital 
to analogue conversion circuitry. 

An object of the invention is 'to provide an 
economical circuit for implementing a DFT. 

Thus we provide an electronic circuit for the 
implementation of a discrete Fourier transform (DFT) , 
which includes input means to which are applied electrical 
signals representing quantities on which the DFT is to be 
performed, and a sum and difference processor (SDP) to 
which those signals are applied . and which can perform 
calculations of the following types ori the signals:- 



. 2 . 0080266 

where x 2 are a sequence of input samples and N is the 
order of the matrix, x is an input signal and U and V 
are the results of the computations. According to the 
invention , in an arrangement as set out above/ the 

5 results generated hy the SDP are collected in an 

addressable buffer store, and those results as assembled 
in the addressable buffer store are applied to a sum of 
products multiplier (SPM) which evaluates expressions 
in the following fqrma.t, quoted for an example in 

10 which N = 5; * ! 
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Yr and Yi are, : 


real and imaginary 


results 



respectively of the DFT, where a^, b^ g i> 2 

are the real and imaginary Fourier coefficients, and 



where 



25 ai + jbl = exp ( jl^) V a~ 2 + j b 2 = exp ( j^) 

Tr 0 = Ur 0 + Ur l + Ur 2 

Ti 0 = Ui 0 + Ui l + Ui 2' 
the resulting output signals from the SPM being 

representative of a DFT of the quantities represented 

by the electrical signals applied to said input means. 

A transmultiplexer is a circuit which acts as 

a bidirectional interface between a set of PCM 

highways and a collection of FDM channels. Thus with 

such a circuit it is also necessary to be able to 

convert from the FDM format into the PCM format, and 

in accordance with the invention this is done by the 
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use of an arrangement which includes filters and an 
inverse DFT ( IDFT ) . The IDFT is implemented in a very 
similar manner to the DFT f except that the values at 
coefficients aire different and _ the sequencing or 
certain data is changed. 

Thus we have a filter arrangement which in 
effect uses the DFT circuit in the reverse sense to 
that used for the PCM to FDM transformation. The 
arrangements to be described herein are intended for 
use in a t ransmultiplexer which interfaces between two 
30-channel PCM systems and five 12-channel FDM groups. 

Embodiments of the invention will now be 
described with reference to the accompanying drawings 
in which: 

Fig. 1 represents schematically a 10 point 
transform factorised into 2 point and 5 point 
transforms. 

Fig, 2 represents schematically a 12 point 
transform factorised into 4 point and 3 point 
transforms. 

Fig. 3 shows a general scheme for 
implementing an efficient Fourier transform. 

Fig. 4 represents data formats used in the 
DFT evaluation. 

Fig. 5 is an outline representation of an SDP. 
Fig. 6 is an addressable buffer store (ABS) 
as used in the arrangement of Fig. 3. 

First we consider the theory of the DFT, and 
show how a number of novel developments lead to ah 
implementation which requires substantially less 
multiplications than a straightforward DFT 
implementation. As such it can be called an FFT. We 
start with a mathematical derivation of the algorithm, 
and follow this with an implementation thereof, 
suitable for implementation using LSI digital circuits. 
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: The Discrete Fourier Transform (DFT) . 
.The .DFT can be expressed as follows :- 
- N^l . P,q,E[l,N-l] 

y p = IZ_x q -wPq ' (i) 

q=0 

where w = or alternatively w = exj> ( j 2 (2) 

("j" denotes the square root of -1). 

Equation, : (1). can be expressed in matrix form, 

as given where N = 5 in equation (3). 

0 0 0 o~o- 
y n w w w w w x n 

0 1 2 3 4 u 
y^ w w w w w x^ . 
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y 0 w. w w w w x 0 

yi 0 3' I 4 2 

y o w w w w w x^ 

ft y| ^ o t 

Y4 w V w w w x 4 
The left-hand side of this is a column vector 
of output samples [y p ]# and the right-hand side is 
the square DFT matrix (referred to as the .W matrix) 
operating on a column vector of input samples 

The elements of the DFT or W matrix can be 
derived as, follows. The value pf_the element in the 
r/th row and £'th column is obtained by evaluating 
expression. (4). 

w (p.q mod N) (.4) 
The power of W/'can be evaluated, modulo N 



since 



W P = W P+N (5) 



In equation (3) and the examples discussed 
the number of points in the data set operated upon 
(i.e. the value of N) is taken to be five. This does 
not imply any restriction on the value of N, and in 
fact different values are used in : the transmultiplexer 
application, referred to. The basic DFT (as given in 
(1) or (3) holds for all values of N. 

„ Equation (3) is the basis for the efficient 

DFT. 

An Efficient Implementation of a DFT 

This technique is applicable to DFT ' s where N 
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is a prime number, but where N is not prime, a known 
technique can be used to factorise the transform into 
prime number transforms, each . implementable as . 
described. The implementation is developed using N=5 
as an example, .although any such prime number, could be 
used. (The case for N«2 is trivial and not appropriate 
to this method) „ , ; - 

The derivation of the efficient transform is 
simplified when N is prime, since all powers of w 
appear once and only once in every row of, the W matrix 
in equation (except for row zero) . This imparts a high 
degree of regularity, to the. subsequent manipulation of 
this equation. 

The terms in the W matrix, (i.e. the powers of 
w in equation (3) can be re-written by noting that 
W P* = w (N-p) p.[l, (N-D/2] (6) 

(w p denotes the complex conjugate of w p ) 
or in particular 

w 4 = w 1 , .w • - w . ; (6b) 
When equation (6b) is substituted in equation 
(3) we get: 

0 0 0 0 0 
y n w w w w w x n 

0 1 2 2* 1* 
y, w w w w w x. 
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X 3 

y^ w v w* w" w** w" 1 " x^ 
Ignoring row 0 of equation (7), all elements in column 
4 of the W matrix in equation (7) are the conjugates of 
elements in the same row in column 1 and likewise for . 
columns 3 and 2 (the leftmost column os column 0). 
This allows some saving in computation* 

The complex multiplications needed. to 
evaluate equation (7) can be taken in pairs (by 
selecting those involving multiplication by a 
coefficient and by its conjugate), and these are now 
examined . 

Let the real and imaginary components of a 
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complex variable be denoted by the upper case 
representation of that variable, with either '_£' or 'i/ 
appended to indicate real or imaginary parts 
respectively. T . ' : : -' 

5 "'-Thus x k = 'Xr k + jXa k ' and c !: =r ; Cr + jci etc. 

Then the sum of a multiplication by a 
coefficient and its conjugate can be expanded as follows: 
- y =• c.x k ; -V- c*V'x^_ k !C ' j ; - 1 • • (8a) 

= : (Cr + : jCi) . (Xr R + ijxi k ) : 
10 +: (Cr' - ^jCi)VX^r N _ k + DXi^k) (8b) 

= ( r dr. (xr k + xr-_ k ) - Ci. (Xi k - Xi N _ R ) ) 
+ j (Cr . (Xi k + Xi^_ k ) •+• Ci. (Xr R '- Xr^_ k ) ) (8c) 



likewise 



z = c*.x k + c . x N _ k (9a) 
15 = (Cr - jCi)' (Xr k + jKi k ) 

+ (Cr + jCi) (Xr N _ k + jXi N _ k ) (9b) 

(Cr. (Xr k - + i xr^_ k ) + Ci . (Xi R -Xi N _ k ) ) • 
+ j (Cr.(Xi k + Xi n _ k ) - Ci.(Xr k -Xr N _ k ) ) (9c) 
Some hew variables can be~introduced to present 
20 the sum and difference terms in equation (8c) and (9c)™. 

u 0 " x 0 

U k " X k + X N-k 

v k = x k " x N-k k e [1 ' (N-D/2] (10a) 
Or expressed in terms of real or imaginary components 
25 Or k = : Xr k + Xr N _ k 

Ui k = xi k + xi N _ k 

Vr k = xr k - Xr N _ k ' (10b) 
vi k = xi k - xi N _ k 
Using these variables, equations (8), (9) can be re-written 
30 as follows :- 

y = [Cr.Ur k - Ci.Vi k ]+ j [Cr.Ui k + Ci.Vr k ] (U) 
z = [Cr.Ur R + Ci.Vi k ]+ j [Cr.Ui k - Ci.Vr k ] (12) 
To make for "a more concise description, a notation is 
introduced for the operation expressed in equation (8), (9). 
35 i.e; y=c # (x k , X N _ k ) = c.x k + c*.x N _ k ' (13) 

and z=c # (x R , x N _ k ) = c*.x k + c.x N _ R (14) 
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The operator can be defined by the following expansion: 

c # (x k , x N _ k ) = [Cr . (Xr k + Xr N _ R ) - Ci.(xi k - xi N _ k )] 

+ j [Cr.(Xi R + xi N -K } + Ci * (Xr k " * r N-k )] (15) 
(i.e;'.as given in. equation (8c) etc. ) 

(The y and z variables are introduced temporarily for 
equations (8) (14), they are used with a different 
meaning in other, equations herein. ) 

Using the results, of equation (8 ) , (9 ) and the 
notation of (13) etc. it is possible to rewrite equation 
(7) as follows :- 

y n w° w° w° U n ,0) 

n i i 
y r w u w 1 w z # (x, ,x 4 ) 

n ' ? l * 

y 2 „ , (x 2 ,x^) (16) 

y 3 

The operator requires four real 

multiplications (as does complex multiplication) , but 
the evaluation of equation (16) required in total 
about half the number of multiplications as does 
equation (7). Thus a significant saving in the number 
of multiplications has been achieved by using the 
relations of equation (6), (9). 

As a refinement, the expressions obtained 
when equation (16). is evaluated can be written down so 
that their terms. are ordered, by the indices of w^. 
Y Q = t 
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Y o 



= x Q + w 1 # (x ir x 4 ) + w 2 # (x 2 ,-x 3 ) 

Y 2 = x Q + v/ 1 *! (x 2 , x 3 ) + w 2 # (x 1 , x 4 ) (17) 

Y 3 = x Q + w 1 # (x 2 , x 3 ) + w 2 # (x 1 

Y 4 = x Q + w 1 # (x lf x 4 ) + w 2 # (x 2 , x 3 ) 

N-l 



where t = / x^ 



q=o 

Finally, for completeness the entire sequence of operations 
necessary to calculate a 5 point DFT is as follows :- 
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First, form the u, v and t (complex valued) terms 



U 0 = X 0 

U l - X l + x 4 



V l = x l r X 4 (18) 
; ' U 2 = X 2 Vx 3 
V 2 = X 2 ~ X 3 

t . = X Q + x x + x 2 + x 3 + x 4 
Let the real and 'imaginary £arts of the powers of w be 
written as : ; : > 

W 2 = a l + jb l 
w = a 2 + jb 2 

Then form the y terms, expressed as real, and imaginary 
components 

Yr o = Tr 0 



Yr x « Ur Q + a 1 .Ur 1 - b i* vi i 

Y^ = Ui Q + "■a- 1 /ni 1 + *> 1 -Vz 1 

Yr 2 = Ur Q + a 1 .Ur 2 + b 1 .Vi ? 

Yi 2 = Ui Q +■ a^;Ui 2 - b x .Vr 2 



4- a 2 .ur 2 - b 2 .Vi 2 

+ a 2 *Ui 2 + b 2 .Vr 2 

+ a 2 .Ur x - b 2 * vi i 

+ a 2 ;ui 1 + b^v^ 



Yr 3 = Dr 0 + a l- Ur 2 " b l ;Vi 2 + a 2 ;Ur l + b 2 - vi i ( 19 ) 
Yi 3 "» Ui Q : + a x .Ui 2 + *> 1 .Vt 2 + a^Uij - b 2 .Vr 1 

Yr 4 = Ur Q + a- 1 .Ur 1 ""+ bj^.Vi^ + a 2 .Ur 2 + b 2 .Vi 2 " 

Yi 4 = Oi 0 + a^U^ - b^V^ + a 2 .Ui 2 - b 2 - Vr 2 

The expressions in equation (19) suggest a 

useful hardware configuration in which a 'Sum of ' 

Products Multiplier' , such as described in our 

co-pending application No. 8132315 (T.J.M. Rossiter 3) 

is loaded with the coefficient values formed by real 

and imaginary components of w , w , w , that is 

1, a^, b lf a 2 , b 2/ and then the data is 

presented to the relevant inputs to calculate the Yr 

and Yip terms. Such a configuration is described 

below. 

The advantages of performing a DFT as described 
are as follows 

1. .. T *?e total number of multiplications needed is 

approximately half that for a normal DFT. 
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2. Re-ordering the data (as in equations (17), 
(19) ) allows very efficient use of a particular form 
of Sum of Products Multiplier, 

3. Performing all the non-trivial multiplications 
in a single step : (i.e. evaluating each line of equation 
(19) completely in a Sum of Products Multiplier) 
produces smaller numerical rounding errors than do some 
alternative FFT algorithms. 

Extending the Efficient DFT to Other Values of N 
We have described how an efficient Fourier 
transform was derived for N points, N being prime. We 
now extend it to the case when N is not prime. The 
techniques used are known and are included to show how 
the ideas can be applied more widely. We do not give a 
general treatment, but instead a worked derivation 
shpwing how the .5 point transform can be extended to a 
10 point transform. 

Where N is not prime it can be factorised into 
a set of factors N Q , N x , n m-i* 

N = Nq x n^ x . . . x N M-1 (20) 

Thus it is possible to factorise the N point 
DFT . into a number of smaller DFTs, each implementable 
using the techniques described above. If the factors 
have no common prime factors they are said to be 
mutually prime, and the resultant factorised DFT has 
particularly attractive properties. But even if the 
factors are not mutually prime, techniques similar to 
those described below can be used, though the count of 
the number of multiplications is higher. 

As an example : N=10 has factors 2, 5 (which 
are mutually prime) . A 10 point Fourier transform can 
be realised as five 2 point transforms followed by two 
5 point transforms, see Fig. 1. A further example is 
given by a 12 point DFT, realised as three 4 point and 
four 3 point transforms, see Fig.. 2. 

The correct sequencing of the input and output 
data is crucial to the operation of the factorised DFT, 
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but since the derivation is complicated, only the 
results are presented here. 

A useful saving in the number of 
multiplications occurs when the transforms are 
factorised. The figures for a regular 10 point DFT,. a 
factorised DFT and at factorised DFT in which the 5 
point factors use the FFT algorithm described above are 
as .follows 

transform No. of Rea l 

Multiplications 
. 10 pt DFT ' ' 304 

Factorised 128 

FFT factors 64 

These figures give the number of non-trivial 
real multiplications i.e. multiplications involving 
factors other than unity. * Thus the FFT algorithm 
described above needs (in this case) less than */4 
the number of multiplications involved in a simple DFT. 

Hardware' implementation of the Efficient DFT 
The logic needed to implement the efficient 
DFT is shown in outline in Fig. 3. The data is first 
processed by a "Sum and Difference Processor" (SDP) 
which can perform calculations such as those of 
equation (10 ) . The SDP "is followed by an "Addressable 
Buffer Store" (ABS) which collects the data generated 
by the SDP and routes it to a "Sum of Products 
Multiplier" (SPM) e.g. as described in our above 
mentioned Application, which evaluates expressions of 
the form of those in equation (19). It can evaluate a 
number of multiplications simultaneously and form their 
sum in one data word period. 

. it is possible to design logic specialised for 

one particular order of DFT (e.g. an SDP and ABS 
designed specifically for a 10 point transform), but 
here we use more general-purpose' logic which is usable 
oh a limited/' tho'ugtV" useful, class of transforms. The 
operations performed by this logic for any particular 
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DFT are then determined by control signals, i.e. the 
logic can be electrically programmed to. suit a number 
of related implementations. ^ 

We first define certain terms arid then 
describe the operation of the SDP and ABS, plus brief 
details of their use in typical DFT computations. 

Definition of Terms, Data Formats etc. 

Abbreviations used are ABS for Addressable 
Buffer Store, SDP for Sum and Difference Processor, and 
SPM for Sum of Products Multiplier. 

The following terms are used herein : 

Part-Word : An individual real data value 
(which may be represented for example by a two's 
complement binary number). 

Complex word : An individual complex data 
value consisting of two Part-Words (representing the 
real and the "imaginary parts) . 

Block : A collection of Complex Words forming 
the input or output data of a DFT i.e. there are N 
Complex Words in the Block of input data for an N point 
DFT. 

Part-Words are composed of M binary bits, M 
being constant but not specified, serially sequenced 
such that the least significant bit occurs first 
followed fcy bits. of . increasing significance. Thus a 
part-wo*£d period is M times the fundamental bit period. 
Two's complement arithmetic is assumed to be used. 
Complex words are normally sequenced with the real part 
preceding the complex part, so that a complex word 
period is twice a part-word period. Finally, blocks 
consist of a ser ially-sequenced collection of complex 
words. The sequencing is specified for a particular 
implementation, but a block period is N times a complex 
word period. These formats are shown in Fig. 4. . 

The fundamental .unit of data in general is the 
complex word, though at times it., may, be necessary to 
manipulate separately the real and imaginary part-word. 
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Likewise the fundamental unit of time is the complex 
word period, though the logic will need synchronisation 
at both the. part-word rate and bit rate. 

Sum and Difference Processor (SDP) 
The SDP shown in outline in Pig. 5 has two 
parallel paths f a' and 'b' which are treated slightly 
differently. The blocks in Pig. 5 have the following 
functions : ^ _ 

The 'Input Gating 1 1 ^can select either^JLri^ut 
data or the stored results of previous computations. 
The selected signal passes to the 'Sum and Difference 1 
circuit 2 which has two outputs. The 'a 1 path output 
is the sum of the two _p^€Es, the 'b 1 path output is the 

difference. . „ _„ 

There is a complex^wbrd^ period delay 3 after 
the sum and difference circuit, which is finnruded- to~ _ 
allow certain data re-ordering. The box marked 'xj* in 
the 1 b* path indicates multiplication by 'j 1 . which can 
be activated or deactivated. Finally the two paths 
pass to the store 4 where they can be held for further 
computation. The^store can also swap over the paths 
(as needed ijr rceX^aTtTlco^ — _ 1 

The output from the SDP is takefr-^f ter the 
delay stage 3 . It includes a Circuit 5 to re-order the 
real "and imaginary part wibrds so tffki-^calculations of 
the form of equation (11), (12) can be performed easily 
by an'SPM. The output can also be passed inter-naily^to 
an accumulator 6 , whose output is available separately;^ 

The operation of these various components of 
the SDP can be modified by control signals, though 
these are not described in detail. The SDP is 
sufficiently versatile to be used in the computation of 
: a number of different DPT 1 s , three examples of which 
are given above. 

Addressable Buffer Store" (ABS) (Fig. 6) 
The purpose of 'the Addressable Buffer Store 
(ABS) is to collect a block of data and then allow a 
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selection of the complex words in that; block to be 
routed to a number, of output lines. The store is so 
configured that one part can be collecting data whilst 
the other part delivers data collected during a 
previous block- At the start of the next block, the 
part which was previously collecting data can be used 
to deliver it, and the other : part can in turn collect 

fresh data . . --- : - ;: ^ 

The ABS, Fig. 6, consists of two main stores 
numbered 0 and 1 and logic to load the stores and route 
the contents of selected locations to the ABS outputs. 
Each store can operate relatively independently and is 
internally partitioned into four banks. The facility 
of loading a store whilst reading previously loaded 
values is realised by using some banks for loading and 
others for reading ... 

The loading of both stores is similar. During 
each complex word period one bank of store 0 and one 
bank of store 1 may be loaded, the banks being , 
independently selected. When a store is loaded the 
existing contents (one complex word per location) •move 
down 1 (e.g. the complex word in 3r moves to 2i and is 
replaced by that in 3i) and the input is loaded into 
the topmost location (e^g. 3i). When a store is not 
selected for loading, the contents remain in the same 
location. Each location holds a complex word. 

There is a slight difference between the two 
stores in the method of reading data. In store 0 one 
location from one bank is selected and fed to output 
Y Q . In store 1 six locations, all from the same 
bank, can be selected simultaneously. The locations 
and outputs are.grouped in three pairs (the "i" and "r n 
locations being grouped together) and the pair of 
locations selected for each output is independently 
controlled. Each "i" output (Yi^; Yi 2 , Yi 3 ) has 
an inversion facility, used to effect the sign 
inversions necessary for complex multiplication. The 
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locations accessed by th.e output selectors can be in 
any of the four banks, the actual^ bank selected being 
deterinijved^by^ an additional contrbl^e<^ninpn tojboth 
stores) . 

By using the vaTiote^ store ac?c r ess 'control^ 
described 1 above, it is possible ^_to ujejftjie— AB4_ to route -'K 
data to_a "sunroof products 1 multiplier, and so compute 
expressions of the f orjn 'given in equat ion Ci9 ) . it is ^ 
also possible touse the ABS^a^a^Time^fo^^^-^^; ' - 
Interchanged", an example of such use' could be in 
arranging the data in^'e ' speci«r^equence required by • '■ 
a particular DFT. - '• '- • - "• ' • 

Examples- of dft'-s ; '- '• . ' .~~ '■ 

The implementation "of three FFT's using the 
SDP, ABS and SPH : devices is described. It is possible — " 
to apply the SDP/ ABS and : £PM to a wider range of DFT's 
than given in these examples, and ^in" fact in the trans- 
multiplexer application, - it was a 14 -point application J 
made up of seven two-point devices feeding two seven 
point devices. In some cases this :: may require a :: 
multiplicity of the three basis devices and/or an ABS 
device with a greater 1 storage ; capacity; : ; • 
■ Example 1. A 5— point DFT ; 1 :t ? 

This is the simplest case : a 5 point - 
transform is performed using the techniques referred to 
above. The SDP evaluates the expressions given 1 in 
equation (18) and generates the u, v and f complex 
terms during a block period. They are then stored in 
the ABS in the following manner :-' 

u 0' t store 00 (or 01 during alternate blocks) 

u l' v l' u 2 ' v 2 store 10 11 during alternate _blocks) 

The SPM is loaded with the <real and imaginary) " ^"" r 

coefficient values for w° f w 1 , „ 2 and then the u ^ 
v and t terms are routed to the SPM to yield the 
expressions given in equation (19). 

Example 2. A 10 point DFT 

in- this case the lO point- transform is factorised 
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into 2 point and 5 point transforms (in the manner shown 
in Fig, 2.1).,. The SDP evaluates the five 2 point DFT 
factors and the sum and difference terms in the 5 point 
factors. The calculations performed by the SDP are given 
in equations (21) , (22) • 



x o 
x o 


+ 
+ 


X 5 
x 5 


Z 0r 
z l> 


<\ 


Z 0 - u 0 
z l " u l 


x 2 

X 2 
X 8 
X 8 


+ 
+ 


X 7 
X 7 
X 3 
X 3 


z 2' * z 2 + . 

2 3' -^ >< ^ Z 2 ~ 

z 9' * 2 3 " 


Z 8 . 
2 8 
Z 9 
Z 9 


- u 2 

- V 2 : 

- u 3 

- v 3 



x 4 


+ 


x 9 


Z 4' 




2 4 




2 6 




U 4 




X 4 




x 9 


: 5 'x 




2 4 




Z 6 




V 4 




X 6 


+ 


X l 


Z 6' 




Z 5 


+ 


Z 7 




U 5 




X 6 




X l 


Z 7' 


-> 


2 5 




z 7 




V 5 


(21) 


fc 0 




*0 


+ z 2 + z 4 


+ 


2 6 




2 8 








fc l 




z l 


+ Z 3 + 2 5 


+ 


z 7 


+ 


2 9 






(22) 



The SDP interleaves calculations for the 2 point 
transforms with the pre-processing for the 5 point 
transforms. The data generated by the SDP is stored in 
the KBS in the following manner 

u Q , t Q store 00 (or 02 during alternate blocks) 

u 1 , t 1 store 01 (or 03 during alternate blocks) 

u 2' v 2' u 4' v 4 store 10 (° r 12 during alternate blocks) 
u 3' v 3' u 5' v 5 store 11 (° r 13 during alternate blocks) 
The five point transforms can then be evaluated in a 
similar manner to the first example (using (u^, u 2 , 

V 2' u 4' v 4 } and (u l' u 3' v 3' u 5' V 5 )J 

in expressions of the form of equation (19). 

Example 3. A 12 point DFT 

This is factorised into three 4 point and four 
3 point transforms, as shown in Figure 2.2. . One SDP 
evaluates the three 4 point DFT factors; the three 
point DFTs can either be evaluated directly by an SPM 
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or a second SDP and SPM with fewer inputs. The 
calculations performed by the first SDP are given in 
equation (23), and those by the second SDP in equation 
(24), (25). 

X 0 + X 6 s 0 : ^ 

X 0 " X 6 S l 
x 3 + x 9 s 2 
X 3 ~ X 9 S 3 




s, - js 3 



10 



x 4 + X 10 

X 4 ~ X 10 
x ? + x x 

x ? - x x 




s 4 + S 6 



- s, 
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20 
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z 8 


X 8 " 


x 2 


s 9 


S 8 


" ; s 10 




z io 


x ll + 


x 5 


S 10^< 
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+ ^ s ll 
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(23) 
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(24) 



(25) 
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The output of the SDP is stored in the.ABS in the 
following manner : - 

u O' u l' U 2" u 3 store OO 

t Q/ t ir t 2 , t 3 . store Ol 

u 4' y 4 ' u 5 / y 5 .. .store 10 

u gf Vg. 7 ' u ? , store 11 

Finally the SPM computes expressions of the form:- 
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* r o 
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Di 4' 




b l* Vt 4 


Yr 8 






+ 


a r 


Ur 4 




b l* Vi 4 


Yi 8 






+ 


a l- 


Ui 4 




b l- Vr 4 



(26) 
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CLAIMS : 

1. An electronic circuit for the implementation 

of a discrete Fourier transform (DFT) , which includes 
input means top which are applied electrical signals 
5 representing quantities on which the DFT is to be 

performed, and a sum and difference processor (SDP) to 
which those signals are applied and which can perform 
calculations of the following types on the signals 

• -• v = x <>--^ •••••• • •• - : • 

10 u k = x R + x N _ k 

V k " X k " X N-k 

where are a sequence of input samples and N is the 
order of the matrix, x is an input signal and U and Vare 
the* results of -the confutations, characterised in this, that the 

15 results generated by the SDP are collected in an 

addressable buffer store, and that" the results as 
assembled in the addressable buffer store are applied 
to a sum of products multiplier (SPM) which evaluates 
expressions in the following format, quoted for an 

20 example in which N = 5; 
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where 


Yr 


and Yi are 


real and imaginary 


results 


respectively of the 


DPT 


, where a 


1' b 


1' 


a 2 


' b 2 



are the real and imaginary Fourier coefficients, and 
where 

35 ax + jbl = exp (j!^)' a 2 + 3 b 2 = exp (j2_) 



Tr 0 = Ur 0 + 0r l + Ur 2 
Ti Q = Ui Q + Ui x + Ui 2 , 
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the resulting output signals from the SPJM being 
representative of a DFT of the quantities represented, 
by the electrical signals applied to said input means. 

2. A circuit arrangement f or .converting the signals 

5 from a number of PGM channels into a frequency division 
; mulitplex (FDM) group, characterised by means to 
convert the PCM words into their linear representations, 
connections over which said, linear representations are 
applied to the said input means, and a polyphase filter 

10 to which the said output signals are' applied. 

3. An electronic circuit for the implementation of 

an inverse Fourier transform (IFT) , which uses a circuit 
as claimed in claim 1, modified in that it operates in 
the reverse direction, i.e. effecting the conversion 
15 from PCM to FDM. 
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